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Abstract 



For the diagonalization of the Hamilton matrix in the Heisenberg model relevant dimen- 
sions are determined depending on the apphcable symmetries. Results aie presented, both, 
by general formulae in closed form and by the respective numbers for a variety of special 
systems. In the case of cyclic symmetry, diagonalizations for Heisenberg spin rings are per- 
'-^ I formed with the use of so-called magnon states. Analytically solvable cases of small spin 

Ch ' rings are singled out and evaluated. 
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5!^ ' 1 Introduction and summary 

-(— > ' 

B- 

' I As far as their magnetic behaviour is concerned, some recently synthesized mole- 

qj ■ cules like "ferric wheels" of six, eight or ten iron ions of spin | [1-4] appear as a 

O . limited array of localized single-particle spins which are adequately described by 

the Heisenberg model [5-8]. The calculation of key quantities like for example the 

spin- spin correlation function becomes easy once the Heisenberg Hamiltonian has 

rS I been diagonalized. For a straight forward diagonalization the dimension d of the 

c3 I Hilbert space H, which for instance isd = {2 s-fl)^ for A^ spins with spin quantum 

number s, may, however, appear prohibitively large even for rather small systems. 

But the obvious symmetries allow to reduce the problem to a set of less ample 

problems according to a decomposition of H into a set of mutually orthogonal 

subspaces. 
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In the case of spin arrays the following symmetry operators may be employed, 
depending on the degree of symmetry exhibited by the system, i.e. the Hamilton 
operator under consideration: 

• the 3 -component of the total spin in the case of axial symmetry, 

• the total spin in the case of full rotational symmetry, 

• the cyclic shift operator for rings of identical single spins with translationally 
invariant coupling, 

• further discrete symmetries like reflection of the spin ring, reflection of the spin 
orientation and reflection of time order. 

Group theoretical arguments are already used to calculate the spectra of some 
molecules by means of an irreducible tensor operator approach [6,9]. In the present 
article an alternative technique is devised. Also, general expressions for the rele- 
vant dimension are derived without taking recourse to any specific procedure of 
diagonalization. 

We will generalize considerations by D. Kouzoudis [10,11] who found analytical 
solutions for special cases of small Heisenberg rings. To this end we are concen- 
trating on the first three symmetries and denote the corresponding simultaneous 
eigenspaces by 7i(S', M, p). After introducing our notation in Sec. 2, we will inves- 
tigate the various symmetries and their corresponding subspaces in Sec. 3. The di- 
mension of the 7i(S', M, p) can be calculated exactly and we derive an explicit for- 
mula for the dimension of the subspaces 7i{S^ M) and 7i(S', M, p), which enables 
the evaluation of the relevant dimension d^. In Sec. 4 we provide recursion formu- 
lae to calculate the Hamilton matrix restricted to the subspaces 7i(S', M, p). To this 
end tools from solid state physics, namely /-magnon states, are utilized [12-14]. 
Section 5 is devoted to examples and applications. For rings of A^ = 5, s = 1 and 
A^ = 8, s = I an exact diagonalization of the Heisenberg Hamiltonian with nearest 
neighbour interaction has been performed using Mathematica® [15]. Further we 
compare the analytical anti-ferromagnetic ground state energies for A^ = 2 ... 8 
with the asymptotic Bethe-Hulthen formula for A^ — > oo [16,17]. 



2 Heisenberg model 



The most general Hamilton operator considered in these investigations consists of 
a spin-spin interaction term Hq and a term Hp describing the interaction with an 
external field which is assumed to point in the 3 -direction 

N N 

h = Ho + Hf = -J2 -Ji^^y) l(^) -K^) -E f-^Bs'{x) . (d 

x,y X 



X, y G G = {1, . . . , A^} label the spin sites modulo A^ and the upper right in- 
dex denotes the spin component. All spins are defined to be dimensionless vector 
observables with the commutation relations 



s''{x),s''{y)\=ieabcS%^)^ccy ■ (2) 

The eigenvalue of (S{x)y is s(x) {s(x) + 1), s{x) being the individually fixed quan- 
tum number at site x. The possible eigenvalues of any s^ (x) are the magnetic quan- 
tum numbers rrix. One may introduce the ladder operators 



s^{x) = s\x)±is^{x) (3) 

to reformulate ^o as 

^o = - E 'Ji^^y) |s'(^)£'(^) + ^[£^(^)S"(^) + £'(^)£^(^)]| • (4) 

x,yeG ^ ^ ' 

The total Hilbert space is spanned by the product basis of the single-particle eigen- 
states of all s^[x) 

s^{x) I mi, . . . , m^, . . . , mjv ) = m^ I mi, . . . , m^, . . . , mjv ) , (5) 

which will be a reference basis in our investigations. 

3 Symmetry operations and relevant dimensions 



Our most general Hamiltonian (1) allows only for invariance with respect to rota- 
tions about the 3-axis. Obviously, Hi^ as well as Rf commute with the 3-component 
S"^ of the total spin 



, §' = E i{^) , (6) 



H,S; 
the possible eigenvalues of which are the total magnetic quantum numbers 



N 



M = -5'n,ax, -•S'max + !,•••, 'S'max with S'^ax = E ^(^) C^) 



x=l 



being the maximum total spin quantum number. The total Hilbert space H for the 
Heisenberg model is the direct sum of all eigenspaces H{M) of S^ 



+'S'max 

7i= n{M). (8) 

M=-5max 

The problem of diagonalizing if in 7i is thus broken up into the corresponding 
problems in each of the 7Y(M). This reduces the dimension one has to cope with 
from 



TV 

d = dim (H) = [] (M^) + 1) (9) 

x=l 

to the respective dim (H{M)). The relevant dimension dr then is the maximum of 
those dimensions, if no other symmetries can be exploited. For given values of M, 
N and of all s{x) the dimension dim (7i(M)) can be determined as the number of 
product states (5), which constitute a basis in 7i(M), with Y.xeG'i^x = M. The 
solution of this combinatorial problem can be given in closed form. 

Theorem: 



dim in{M)) 



(^max-M)! 




Smax-M N i _ 2s(x)+l 

n / . 



x=l 



(10) 



z=0 



For equal single-spin quantum numbers s(l) = ■ ■ ■ = s{N) = s, and thus a maxi- 
mum total spin quantum number of S'max = Ns, Eq. (10) simplifies to 



dim{niM)) = f{N,2s + l,S^^^-M) with (11) 

n=0 \' 



n \ N - 1 



In both formulae (10) and (11), M may be replaced by \M\ since the dimension of 
l-i{M) equals those of l-i{—M). [f/ jJi] in the sum symbolizes the greatest integer 
less or equal to f/ fi. Eq. (1 1) is known as a result of De Moivre [18]. 

Proof: The proof of (10) and (11) may be accomplished by comparing any product 
state (5) with the completely aligned state 

I f] ) = I mi = s(l), ma = s(2), . . . , m^ = s{N) ) , (13) 

which is also called magnon vacuum state, see next section. Evidently it is also an 
eigenstate of H and in the ferromagnetic case a ground state of Hq. Using the pos- 
sible decrements a(x) = s(x) — m^. of the magnetic quantum number at each spin 
site we want to construct a generating function F(z) of the numbers dim (7i(M)) 
as a polynomial in z 



25max 

F{z)= J2 dim {n{S^,.-n))z^. (14) 

n=0 

The consideration, that the number of product states with fixed magnetic quantum 
number M corresponds to the related number of sequences of decrements, leads to 
the following definition 



\ ^ 1 _ z^s{x)+l 

^w = n E ^"^^M = n ._, • (15) 

x=l \a(x)=0 j x=\ i ^ 

Eq. (10) is then obvious and (11) is easily inferred from (10) by using Leibniz's 
theorem [19] to calculate multiple derivatives of the product {\—z^^^^)^-(\—z)^^. 



In the special case of identical single-particle spins, Hp is proportional to S^ and 
therefore commutes with i/o- The simultaneous eigenstates of Hi^ and S^ are then 
also eigenstates of Hp. Hp causes a splitting proportional to M of the eigenvalues 
of Hq, which otherwise would not depend on M. 

In what follows, consideration is alternatively restricted to one of two special cases. 
In the first case (I) the Hamilton operator is field free, H = Hq, and thus invariant 
under any rotation. In the second case (II) all individual spins are identical. 

The total spin quantum number will be denoted by S. It has values in the set 

{0 < 5'min, 5'min + 1, • • • , 5'max} and 5^ has eigenvalues S{S + 1). The simulta- 
neous eigenspaces 7Y(S', M) of S*^ and S^ are spanned by eigenvectors of H. The 
one-dimensional subspace 7i(M = S'max) = '^(•S'max, •S'max), especially, is spanned 
by I f2 ) . The total ladder operators are 

S^ = S^±iS^. (16) 

For S > M, S~ maps any normalized if-eigenstate e H{S,M + I) onto an H- 



eigenstate G n{S, M) with norm ^J S{S + 1) - M{M + 1). 

Lemma: For < M < S'max, 7i(M) can be decomposed into orthogonal subspaces 

n{M) = n{M,M)®S'n{M + i) (iv) 

with 



s-n{M + i)= nis,M). (18) 

S>M+l 



Proof: The 7i(S', M) represent an orthogonal decomposition of Ti.{M). It is there- 
fore sufficient to identify the orthogonal complement of S~'H{M+1) within Ti^M) 
as n{M,M). This, in turn, is clear if any \i)) ^ G n{M) with | V^ ) ± 
S~H{M + 1) vanishes on application of S^. The defining property of all such 

I ■?/') is 

(^lg-|0) = o y\(t)) en{M + i). (19) 

But then S+\^) e n{M + 1) and ( S+tp | ) = which proves g+ | ^ ) = 0. ■ 

In consequence, the diagonalization of H inH has now been traced back to diago- 
nalization in the subspaces H{S, S), the dimension of which are for S < S^ax 

dim {n{S, S)) = dim {n{M = S)) - dim {n{M = S + 1)) (20) 

and can be calculated according to Eq. (10) in case (I) or Eq. (11) in case (II). The 
relevant dimension dr is the maximum of the numbers given by Eq. (20). 

For many of the systems of interest here, the individual spin quantum numbers are 
equal, the coupling coefficients J have the special property J{x,y) = J(\x—y\) and 
the applied magnetic field is homogeneous. This renders the Hamiltonian invariant 
with respect to translations on the group G = {!,..., A^} of spin sites . Any such 
translation is represented by the cyclic shift operator T or a multiple repetition. T 
is defined by its action on the product basis (5) 

T I mi, ... , rriN-i, rriN) = \ rriN, mi, ... , rriN^i ) ■ (21) 

The eigenvalues of T are the A^-th roots of unity, z = exp{— ip}, where p will be 
called magnon momentum and can take the following N values modulo 2n from 



G = L\p=^,k = 0,...,N-l\ , (22) 

where k will be called translational quantum number. Clearly, T commutes now 
with both Hamiltonian and total spin. Any H{S, M) is decomposed into simultane- 
ous eigenspaces 7i(S', M,p) of S*^, S^ and T, and diagonalization can be reduced 
to a diagonalization in the subspaces H(S', S^p). The according reduction of the 
relevant dimension can be quantified if the values of p and their degeneracy within 
?i(S', S) are known. As a rule of thumb one finds a tendency for equal degeneracy 
and thus 



dr^^ max dim {n{S, S)) . (23) 



It is, however, possible to determine rf,. exactly for any A^ and s. To this end we 
introduce the notations of cycles, proper cycles and epicycles. A special decompo- 
sition of H into orthogonal subspaces can be achieved by starting with the product 
basis (5) and considering the equivalence relation 

1^)^ 10)^ |^)=T"|0),nG{l,2,...,iV} (24) 

for any pair of states belonging to the product basis. The equivalence relation then 
induces a complete decomposition of the basis into disjoint subsets, i.e. the equiv- 
alence classes. A "cycle" is defined as the linear span of such an equivalence class 
of basis vectors. The obviously orthogonal decomposition ofH into cycles is com- 
patible with the decomposition ofH into the various Ti.{M) but not, generally, with 
the decomposition of 7i(M) into its subspaces 'H{S, M). Evidently, the dimension 
of a cycle can never exceed N . Cycles are called "proper cycles" if their dimension 
equals N, they are termed "epicycles" else. One of the N primary basis states of a 
proper cycle may arbitrarily be denoted as 

I ^1 ) = I mi, . . . , m^;, . . . , mAT ) (25) 

and the remaining ones may be enumerated as 

|^„+i)=T"|^i),n = l,2,...,iV-l. (26) 

The cycle under consideration is likewise spanned by the states 

»> = ^=E (e-'-^T) |«v+i> (27) 



;^E(e-*l)I^V« 



which are eigenstates of T with the respective magnon momentum p( A;) = 2Txk/N. 
Consequently, every k (every p{k)) occurs once in a proper cycle. 

The reader will easily verify the validity of the following related statements: 

• An epicycle of dimension D is spanned by D eigenstates of T with each of the 
translational quantum numbers /c = 0, N/D, . . . ,(D — 1)N/D occurring exactly 
once. 

• For a primary product state in an epicycle of dimension D, the spin sites are 
grouped into D "subrings", each subring having a constant magnetic decrement 
s — rrix =: a on all of its places. The subring forms a pattern which repeats itself 
on a cyclic shift of D steps. The total number of spin sites in a subring is N/D. 

• Therefore, for M < S'max, 7i(M) contains a cycle (epi or proper) of dimension 
D if and only if D is a divisor of N, including D = \ and D = N, and D^S^ax — 
M)/N is an integer. The set of such allowed D will be denoted as C{N, s, M). 



Let K.m{D) be the linear span of all cycles of dimension D occurring within H{M) 
and let ks{N, M, D) denote its dimension. According to the above remarks we have 



dim{l-i{M))= J2 ks{N,M,D) (28) 

DeC{N,s,M) 

= ks{N,M,N)+ Y. ks{N,M,D). 

DeC{N,s,M);DytN 

Because to each epicycle of dimension D there corresponds exactly one proper 
cycle with D spin sites and total magnetic quantum number MD/N, and vice versa, 
(28) may be rewritten as 



dim{n{M)) = ks{N,M,N)+ Yl ks{D,^,D), (29) 

DeC{N,s,M) 

which, together with (11), may be used as a recursion relation for the function 
ks(N, M, N) and hence also for ks(N, M, D). This recursion relation can be trans- 
formed into an explicit formula which reads 



J(N.s.M) ^ ^ \ 



h{N,M,N)= Y ^(n)/lA2s + l,^^%^) (30) 

D€C{N,s,M) ^^^ ^ ^^ 



where / is taken over from Eq. (12) and 



, , ,,-1)'" if z/ is a product of m different primes, ^,.^ 

"'"'^^O else. <"' 

The proof rests on some elementary combinatorics and will be omitted here. Further 
we have for the simultaneous eigenspaces 7i(-, M, p{k)) of S*'^ and T 



dim(7^(-,M,p(A;)))= f] A; JiV, M, ^) ^ (32) 

u=l 

u\N;u\k 

and 



Aim{n{S,M,p)) = Aim{n{-,S,p))-dim{n{-,S+l,p)) (33) 

for < S* < S'max and \M\ < S. This allows the explicit calculation of the relevant 
dimensions for any given A^ and s, for examples see table 2. 



4 Diagonalization of the Hamiltonian in the magnon basis representation 



In this section we consider the Hamiltonian (1) with J(x, y) = J{\x ~ y\), B = 
and equal single spin quantum numbers 

H = - J2 Ji\x-y\)six)-siy), yxeG:six) = s. (34) 

In order to calculate the matrix elements of this Hamiltonian restricted to subspaces 
Hi^S, M, p) it is recommendable to use a basis of vectors which are already adapted 
to the problem. We found it most convenient to work with a basis constructed from 
so-called magnon states, used for example in [14]. These magnon states should 
not be confused with those defined by Holstein and Primakoff [13]. The pertinent 
definitions are the following. 

For any function / of spin sites x E G the discrete Fourier transform is as usual 
defined by 



/p ^= ^ E e^'' /(^) ' P e G . (35) 

The analogous transformation may be applied to linear operators A{x), where A{x) 
for instance may be a single-particle spin component at site x 



§^=^E e^"^f(^) (36) 

or a ladder operator (3), which yields 

Then the Hamiltonian (34) can be written in terms of the S^ and §p as (see [14]) 
B = -^Y.JA-Iv (38) 

pad ^eG 

It is straight forward to obtain the following commutation relations 



c+ q+ 



o- q- 
x,p ' lea 



0, 



iV' 



S[ 



p+q 5 



Co C± 



1 



iV' 



5, 



p+f? 



(39) 
(40) 



from the corresponding commutation relations of the spin operators at site x. 
We define "/-magnon states" with momenta pi, . . . ,p; as 



fip) 



n 



pi,---,pi 



) •" X.P1 §Pi ' ' ' ^Pi I ^^ ) 



(41) 



for pi, E G,v = 1, . . . , / and / = 0, . . . , 2Ns, \ fi ) being the magnon vacuum 
state as given in Eq. (13). For / > 2Ns \ fi^ ) will be the zero vector. Since the Sp^ 
commute, a unique representation of (41) may be achieved by postulating pi < 

P2<---<Pl- 

One-magnon states | fip ) are orthogonal 



\\Lp \ \Lq / ZSOpq , 

which follows from (40) and 



(42) 



Si \n) = ^s6pfl \^) . 
Moreover, the one-magnon states are already eigenvectors of H and T 



(43) 



H\Vtp)=Ep\Vtp) , Ep = 2s Jp + s{Ns - 2) Jq 



T\np) 



-ip 



fip). 



The last relation may be derived from 



T^ S^ = e-^P^ S^ T^ Vx G G, Vp G G . 



(44) 

(45) 



(46) 



For / > 1 the /-magnon states are in general no longer orthogonal or linearly in- 
dependent. But they span the total Hilbert space H. and they behave nicely with 
respect to the ladder operator S^ , namely 



and they span the subspaces H^M, p), which follows from 



(47) 



10 



T I fipi,...,p, )=ey.^l-i^Pv\ \ ^pj,...,pi ) (48) 

Si I n,,_,, ) = iNs- I) I n,,_,, ) . (49) 

Let the magnon states | f]^ ) be ordered as follows 

I fi ) (50) 

I ^0 )? I ^1 )) • • • ? I ^Af-l ) 



where each row is ordered lexicographically. We will call the basis of H generated 
by applying the well-known Gram-Schmidt orthonormalization procedure to the 
sequence (50) "magnon basis". Note that the Gram-Schmidt procedure will some- 
times produce zero vectors, because the | Clf ) are linearly dependent. The main 
result of this section is then the following proposition. 

Proposition: The magnon basis contains a subbasis for each subspace H^S, M, p), 

S = Ns,Ns-l,...,\M\; peG. 

Proof: Since the | fi^) span the subspaces l-i{M,p) and different l-i{M,p) are 
orthogonal, the magnon basis spans 'H{M,p), too. The subspace with maximal 
S = Ns is spanned by the magnon states 

\fi), I ^0 )) I ^00 ); • • • I ^(2Afs)zeros ) ; (51) 

which also occur in the magnon basis up to normalization. These vectors also span 
the one-dimensional spaces 7^(5*, M,p) with S = Ns; M = —S, . . . , S;p = 0. 
Hence the proposition holds for S = Ns. We then may proceed by induction. 
Under the assumption that the proposition holds for all S > Sq, we have to show 
that it holds for S = Sq. Consider the subspaces l-i{M) with M = Sq. According 
to (18) we have 

n{So) = s-niSo + i)(BniSo,So) . (52) 

Hence 7i(5'o, Sq) is spanned by vectors from the magnon basis and the same holds 
for every H{So,So,p),p E G. For the other subspaces H{So,M,p), M < Sq, we 
have 

n{So,M,p) = {s-y"'"^ niSo,So,p) (53) 

and we know that the restrictions 



11 



S- : n{M + 1) — ^ n{M) , M > 0, (54) 



are isometrics up to a factor JS{S + 1) — M{M + 1). Hence S^ maps vectors 
from the magnon basis onto vectors from the magnon basis (up to a factor), which 
concludes the proof. ■ 

We are left with two tasks. First, in order to calculate the magnon basis, we need a 
formula for the inner product between magnon states. We have no explicit formula, 
but using the commutation relations (39) and (40) we derived a recursion relation, 
which is sufficient for computer algebraic calculations. Similarly, we can express 
the matrix elements of the Hamiltonian between two magnon states by finite sums 
containing only inner products of magnon states, which completes the second task. 
The results are 



v=l 

2 
~T7 Z^ \ ^^P2v,Pi I ^^<?lv■■,l3I'v■■,(<?M^^(9M+'^!^-Pl))v■■,g^ / ' 

where the symbol q^ denotes deletion of the index q^ and (g^ ^^ {Qf^ + Qu — Pi)) 
denotes replacement of the index g^ by the given expression which has to be under- 
stood modulo N. Further 



\ "pi,...,p; I -H I iigi,...,g; / 
~ 2^ "^P \ ^^P,Pl, — ,Pl I ^^P.9lv,'?i ) ' "^0 { "pi,...,p; I "<Ji,...,IJ; / I -5 ^V — ZSl 

ped ^ 

~ 1 ' 

+ Z^ "^p "/Y Z^ \ ^^piv,(gi'^('?i'+p)),---,Pi I ^^gi,---,(gM^('?M+p))v,'?i / • 

Now we have all what is required in order to write a program [15] which reduces the 
diagonalization of H to the diagonalization of the submatrices of H with respect to 
the magnon basis of the relevant subspaces 7i(S', M,p). The final diagonalization 
can be done numerically or in the few cases, where the relevant dimension is less 
or equal 4, analytically. Table 2 shows in bold face for which cases analytical solu- 
tion will be possible. Two examples for eigenvalues and multiplicities are given in 
the next section. The complete routine, which also presents the eigenvectors, may 
be downloaded from the internet [15]. It may be possible to exactly solve other 
cases if the discrete symmetries mentioned in the introduction are taken into ac- 
count. The cases A^ = 2, 3, 4 can, in principle, also be solved by Clebsch-Gordan 
decomposition of the Hilbert space. 
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Table 1 

Dimensions of the subspaces 7^(5, M,p(fc)) for A'^ = 6, s = ^. Each column can be created 
by applying the ladder operator (16) yielding {2S+ 1) -dimensional subspaces. One realizes 
a discrete symmetry between columns for k and 6 — k. 

5 Examples 



5.1 Relevant dimension of spin rings 



In order to illustrate the above considerations let us discuss the example A^ = 6, 
s = |, which has been solved by D. Kouzoudis [11]. The decomposition (8) yields 
the following relation between the dimensions of the involved Hilbert spaces 



dim (?i) = 2^ = 1 + 6 + 15 + 20 + 15 + 6 + 1 = ^ dim (?i(M)) .(57) 

M=-S'max 

Obviously, | ^ ) = | + + + + ++ ) spans the one-dimensional space 7i(M) with 
M = S'max = Ns = 3, which is an epicycle of dimension D = 1. l-i{M = 2) 
consists of the proper cycle generated e.g. by | — h + + ++ ) and hence is six- 
dimensional. It is spanned by the six eigenstates of T and H given in Eq. (27), 
which are one-magnon states. The next space 7i(M = 1) consists of two proper cy- 
cles generated e.g. by | h + ++ ) and | — I h ++ ), respectively, and one 

three-dimensional epicycle generated by | — \- -\ 1-+ ). Hence dim (7i(M = 1)) 

6 + 6 + 3 = 15. Note that the last mentioned epicycle corresponds to the proper 
cycle generated by | — h + ) with N = 3, M = ^. The largest space 7i(M = 0) 

is spanned by three proper cycles generated by | 1 1-+), | + H 1 ), 

I + + H ) and one two-dimensional epicycle generated by | H 1 1 — ). 
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24 
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120 

24 


609 
105 


3150 
450 
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88900 
9884 


484155 
48483 



Table 2 

Relevant dimension assuming only invariance with respect to rotations (upper rows) and 

assuming also invariance with respect to cyclic shifts (lower rows). 

Consequently, dim(7i(M = 0)) = 6 + 6 + 6 + 2 = 20. The remaining spaces 
7i(M) with M < have the same cycle structure as 7i(|M|) due to the symmetry 



The dimensions of the subspaces 'H(S', M,p{k)) are given in table 1. Note that in 
accordance with the rule of thumb (23) the numbers dim (H{-, M,p)) are almost 
uniformly distributed with respect to p. The deviations from the uniform distribu- 
tion for \M\ = 0, 1 can be explained by the extra eigenvalues of T produced by the 
mentioned epicycles of dimensions three and two. 

The next table 2 shows the relevant dimensions for a wide variety of N and s. The 
ratio of upper-row to lower-row entries is roughly N in accordance with the rule 
of thumb, Eq. (23). Dimensions less than five are given in bold, these cases can be 
solved analytically. Two examples are presented in the next subsection. 



5.2 Exactly solvable systems 



Analytical solutions are for instance possible in the cases A^ = 5, s = 1 and A^ = 
8, s = i, which to our knowledge have not yet been published. Table 3 contains 
the exact eigenvalues and multiplicities for the Heisenberg Hamiltonian Eq. (34) 
with nearest neighbour interaction and A^ = 5, s = 1 and Fig. 1 shows the same 
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|(6± [V5 + 2V23cos(a)]) 




deg(E) 


14 
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10 


6 




E/J 


1 (6 ± [^/5 - ^/23cos(q)] ± ^/69sm(a)) 






deg(E) 


6 





Table 3 

Eigenvalues and degeneracies for the Heisenberg Hamiltonian Eq. (34) with nearest neigh- 



bour interaction and iV = 5, s = 1. The value of a is a = ^arctan 
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Fig. 1. Eigenvalues and degeneracies for the Heisenberg Hamiltonian Eq. (34) with 
N = 5,s = 1 (l.h.s.) and N = 8,s = \ (r.h.s.). The lines are drawn as a guide for 
the eye. 

information graphically, but also for A^ = 8, s = |. The analytical, but rather 
lengthy expressions for the latter case may be evaluated using a Mathematica® 
script we provide [15]. Here we would like to present only the eigenvalue which in 
the antiferromagnetic case corresponds to the ground state energy 



Eq/J = -<^2 + V13COS 



1 , CiV^" 

-arctan 

3 Is, 



for iV = 8, s 



(58) 



In order to clarify the structure of the energy spectrum we calculate the smoothened 
density of states. To this end the degeneracy of an energy eigenvalue En is first 
divided by the mean distance to the neighbouring eigenvalues, then the resulting 
function is linearly interpolated and convoluted with a characteristic function of 
width 1. This density of states p{E), which is displayed for the Heisenberg Hamil- 
tonian Eq. (34) in Fig. 2 for A^ = 5, s = 1 (l.h.s.) and iV = 8, s = | (r.h.s.), shows 
an interesting band structure. 
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Fig. 2. Smoothened density of states p{E) for the Heisenberg hamiltonian Eq. (34) with 
N = 5,s = l (l.h.s.) and iV = 8, s = i (r.h.s.). 
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Fig. 3. Comparison of analytical {N < 8) and numerical {N > 8) ground state energies 
(symbols) for antiferromagnetic coupling with the large A^ limit (solid line) of Bethe and 
Hulthen [16,17]. 

5.3 Bethe-Hulthen 



The calculated eigenvalues, Table 3 and Figure 1, correspond to the ferromagnetic 
case for J > and to the antiferromagnetic case for J < 0. Thus our data also 
contain the antiferromagnetic ground state energy Eq{N) and we can compare 
Eq{N)/{NJ) with the exact limit 



lim 

7V->oc 



E,{N) 

NJ 



2In(2)-- 



(59) 



known as the Bethe-Hulthen formula [12,16,17]. In Fig. 3 analytical {N < 8) and 
numerical {N > 8) ground state energies are displayed by symbols and the large 
A^ limit by a solid line. It turns out that the limit is approached from above for even 
N and from below for odd A^ and that the approach is faster for even values of A^. 
This shows the effect of spin frustration. 
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